Numerical solution and convergence analysis of steam injection in heavy oil reservoirs     

H. Hajinezhad  Ali R. Soheili  Mohammad R. Rasaei  F. Toutounian 《Computational Geosciences》2018,22(6):1433-1444

In this paper, the numerical methods for solving the problem of steam injection in the heavy oil reservoirs are presented. We consider a 3-dimensional model of 3-phase flow, oil, water, and steam, with the effect of 3-phase relative permeability. Interphase mass transfer of water and steam is considered; oil is assumed nonvolatile. We apply the simultaneous solution approach to solve the corresponding nonlinear discretized partial differential equation in the fully implicit form. The convergence of finite difference scheme is proved by the Rosinger theorem. The heuristic Jacobian-Free-Newton-Krylov (HJFNK) method is proposed for solving the system of algebraic equations. The result of this proposed numerical method is well compared with some experimental results. Our numerical results show that the first iteration of the full approximation scheme (FAS) provides a good initial guess for the Newton method. Therefore, we propose a new hybrid-FAS-HJFNK method while there is no steam in the reservoir. The numerical results show that the hybrid-FAS-HJFNK method converges faster than the HJFNK method.  相似文献   

Sensitivity of the Gravity Recovery and Climate Experiment (GRACE) to the complexity of aquifer systems for monitoring of groundwater     

Yashwant B. Katpatal  C. Rishma  Chandan K. Singh 《Hydrogeology Journal》2018,26(3):933-943

The Gravity Recovery and Climate Experiment (GRACE) satellite mission is aimed at assessment of groundwater storage under different terrestrial conditions. The main objective of the presented study is to highlight the significance of aquifer complexity to improve the performance of GRACE in monitoring groundwater. Vidarbha region of Maharashtra, central India, was selected as the study area for analysis, since the region comprises a simple aquifer system in the western region and a complex aquifer system in the eastern region. Groundwater-level-trend analyses of the different aquifer systems and spatial and temporal variation of the terrestrial water storage anomaly were studied to understand the groundwater scenario. GRACE and its field application involve selecting four pixels from the GRACE output with different aquifer systems, where each GRACE pixel encompasses 50–90 monitoring wells. Groundwater storage anomalies (GWSA) are derived for each pixel for the period 2002 to 2015 using the Release 05 (RL05) monthly GRACE gravity models and the Global Land Data Assimilation System (GLDAS) land-surface models (GWSA_GRACE) as well as the actual field data (GWSA_Actual). Correlation analysis between GWSA_GRACE and GWSA_Actual was performed using linear regression. The Pearson and Spearman methods show that the performance of GRACE is good in the region with simple aquifers; however, performance is poorer in the region with multiple aquifer systems. The study highlights the importance of incorporating the sensitivity of GRACE in estimation of groundwater storage in complex aquifer systems in future studies.  相似文献   

Hydrogeology and sustainable future groundwater abstraction from the Agua Verde aquifer in the Atacama Desert,northern Chile   总被引：1，自引：0，他引：1  

Javier Urrutia  Jorge Jódar  Agustín Medina  Christian Herrera  Guillermo Chong  Harry Urqueta  José A. Luque 《Hydrogeology Journal》2018,26(6):1989-2007

The hyper-arid conditions prevailing in Agua Verde aquifer in northern Chile make this system the most important water source for nearby towns and mining industries. Due to the growing demand for water in this region, recharge is investigated along with the impact of intense pumping activity in this aquifer. A conceptual model of the hydrogeological system is developed and implemented into a two-dimensional groundwater-flow numerical model. To assess the impact of climate change and groundwater extraction, several scenarios are simulated considering variations in both aquifer recharge and withdrawals. The estimated average groundwater lateral recharge from Precordillera (pre-mountain range) is about 4,482 m³/day. The scenarios that consider an increase of water withdrawal show a non-sustainable groundwater consumption leading to an over-exploitation of the resource, because the outflows surpasses inflows, causing storage depletion. The greater the depletion, the larger the impact of recharge reduction caused by the considered future climate change. This result indicates that the combined effects of such factors may have a severe impact on groundwater availability as found in other groundwater-dependent regions located in arid environments. Furthermore, the scenarios that consider a reduction of the extraction flow rate show that it may be possible to partially alleviate the damage already caused to the aquifer by the continuous extractions since 1974, and it can partially counteract climate change impacts on future groundwater availability caused by a decrease in precipitation (and so in recharge), if the desalination plant in Taltal increases its capacity.  相似文献   

Groundwater sustainability and groundwater/surface-water interaction in arid Dunhuang Basin,northwest China   总被引：1，自引：0，他引：1  

Jingjing Lin  Rui Ma  Yalu Hu  Ziyong Sun  Yanxin Wang  Colin P. R. McCarter 《Hydrogeology Journal》2018,26(5):1559-1572

The Dunhuang Basin, a typical inland basin in northwestern China, suffers a net loss of groundwater and the occasional disappearance of the Crescent Lake. Within this region, the groundwater/surface-water interactions are important for the sustainability of the groundwater resources. A three-dimensional transient groundwater flow model was established and calibrated using MODFLOW 2000, which was used to predict changes to these interactions once a water diversion project is completed. The simulated results indicate that introducing water from outside of the basin into the Shule and Danghe rivers could reverse the negative groundwater balance in the Basin. River-water/groundwater interactions control the groundwater hydrology, where river leakage to the groundwater in the Basin will increase from 3,114?×?10⁴ m³/year in 2017 to 11,875?×?10⁴ m³/year in 2021, and to 17,039?×?10⁴ m³/year in 2036. In comparison, groundwater discharge to the rivers will decrease from 3277?×?10⁴ m³/year in 2017 to 1857?×?10⁴ m³/year in 2021, and to 510?×?10⁴ m³/year by 2036; thus, the hydrology will switch from groundwater discharge to groundwater recharge after implementing the water diversion project. The simulation indicates that the increased net river infiltration due to the water diversion project will raise the water table and then effectively increasing the water level of the Crescent Lake, as the lake level is contiguous with the water table. However, the regional phreatic evaporation will be enhanced, which may intensify soil salinization in the Dunhuang Basin. These results can guide the water allocation scheme for the water diversion project to alleviate groundwater depletion and mitigate geo-environmental problem.  相似文献   

Upscaling permeability for three-dimensional fractured porous rocks with the multiple boundary method     

Tao Chen  Christoph Clauser  Gabriele Marquart  Karen Willbrand  Thomas Hiller 《Hydrogeology Journal》2018,26(6):1903-1916

Upscaling permeability of grid blocks is crucial for groundwater models. A novel upscaling method for three-dimensional fractured porous rocks is presented. The objective of the study was to compare this method with the commonly used Oda upscaling method and the volume averaging method. First, the multiple boundary method and its computational framework were defined for three-dimensional stochastic fracture networks. Then, the different upscaling methods were compared for a set of rotated fractures, for tortuous fractures, and for two discrete fracture networks. The results computed by the multiple boundary method are comparable with those of the other two methods and fit best the analytical solution for a set of rotated fractures. The errors in flow rate of the equivalent fracture model decrease when using the multiple boundary method. Furthermore, the errors of the equivalent fracture models increase from well-connected fracture networks to poorly connected ones. Finally, the diagonal components of the equivalent permeability tensors tend to follow a normal or log-normal distribution for the well-connected fracture network model with infinite fracture size. By contrast, they exhibit a power-law distribution for the poorly connected fracture network with multiple scale fractures. The study demonstrates the accuracy and the flexibility of the multiple boundary upscaling concept. This makes it attractive for being incorporated into any existing flow-based upscaling procedures, which helps in reducing the uncertainty of groundwater models.  相似文献   

Tracing submarine groundwater discharge flux in Tolo Harbor,Hong Kong (China)     

Yi Liu  Jiu Jimmy Jiao  Ho Kwan Cheng 《Hydrogeology Journal》2018,26(6):1857-1873

Submarine groundwater discharge (SGD) is an important pathway for groundwater and associated chemicals to discharge to the sea. Groundwater levels monitored along a transect perpendicular to the shoreline are used to calculate SGD flux from the nearshore aquifer to Tolo Harbor, Hong Kong (China). The calculated SGD flux—recharge/discharge measured with Darcy’s Law methods—agrees well with estimates based on geo-tracer techniques and seepage meter in Tolo Harbor during previous studies. The estimated freshwater SGD is 1.69–2.0 m²/d at the study site and 0.3?±?0.04 cm/d for the whole of Tolo Harbor, which is comparable to the river discharge (0.25?±?0.07 cm/d) and precipitation (0.45?±?0.15 cm/d). The tide-driven SGD in the intertidal zone is 13.98–17.59 m²/d at the study site and 2.42?±?0.56 cm/d for the whole of Tolo Harbor. The SGD occurring in the subtidal zone and the bottom of Tolo Harbor is 3.12?±?4.63 cm/d. Fresh SGD accounts for ~5% of the total SGD, while the rest (~95%) is contributed by saline SGD driven by various forces. About 96% of the tide-driven SGD in the intertidal zone occurs in the ebbing tide period because the head difference between the groundwater level and sea level is great during this period. Tide-driven SGD in the spring tide is ~1.2 times that during neap tide. The tidal fluctuation amplitude and tide-driven SGD in the intertidal zone are positively correlated to each other; thus, a spring neap variation of the tide-driven SGD is observed.  相似文献   

Characterising the vertical separation of shale-gas source rocks and aquifers across England and Wales (UK)     

Sian E. Loveless  John P. Bloomfield  Robert S. Ward  Alwyn J. Hart  Ian R. Davey  Melinda A. Lewis 《Hydrogeology Journal》2018,26(6):1975-1987

Shale gas is considered by many to have the potential to provide the UK with greater energy security, economic growth and jobs. However, development of a shale gas industry is highly contentious due to environmental concerns including the risk of groundwater pollution. Evidence suggests that the vertical separation between exploited shale units and aquifers is an important factor in the risk to groundwater from shale gas exploitation. A methodology is presented to assess the vertical separation between different pairs of aquifers and shales that are present across England and Wales. The application of the method is then demonstrated for two of these pairs—the Cretaceous Chalk Group aquifer and the Upper Jurassic Kimmeridge Clay Formation, and the Triassic sandstone aquifer and the Carboniferous Bowland Shale Formation. Challenges in defining what might be considered criteria for ‘safe separation’ between a shale gas formation and an overlying aquifer are discussed, in particular with respect to uncertainties in geological properties, aquifer extents and determination of socially acceptable risk levels. Modelled vertical separations suggest that the risk of aquifer contamination from shale exploration will vary greatly between shale–aquifer pairs and between regions and this will need to be considered carefully as part of the risk assessment and management for any shale gas development.  相似文献   

Impact of model complexity and multi-scale data integration on the estimation of hydrogeological parameters in a dual-porosity aquifer     

Elena Tamayo-Mas  Marco Bianchi  Majdi Mansour 《Hydrogeology Journal》2018,26(6):1917-1933

This study investigates the impact of model complexity and multi-scale prior hydrogeological data on the interpretation of pumping test data in a dual-porosity aquifer (the Chalk aquifer in England, UK). In order to characterize the hydrogeological properties, different approaches ranging from a traditional analytical solution (Theis approach) to more sophisticated numerical models with automatically calibrated input parameters are applied. Comparisons of results from the different approaches show that neither traditional analytical solutions nor a numerical model assuming a homogenous and isotropic aquifer can adequately explain the observed drawdowns. A better reproduction of the observed drawdowns in all seven monitoring locations is instead achieved when medium and local-scale prior information about the vertical hydraulic conductivity (K) distribution is used to constrain the model calibration process. In particular, the integration of medium-scale vertical K variations based on flowmeter measurements lead to an improvement in the goodness-of-fit of the simulated drawdowns of about 30%. Further improvements (up to 70%) were observed when a simple upscaling approach was used to integrate small-scale K data to constrain the automatic calibration process of the numerical model. Although the analysis focuses on a specific case study, these results provide insights about the representativeness of the estimates of hydrogeological properties based on different interpretations of pumping test data, and promote the integration of multi-scale data for the characterization of heterogeneous aquifers in complex hydrogeological settings.  相似文献   

Spatial resolution of transport parameters in a subtropical karst conduit system during dry and wet seasons     

Anna Ender  Nadine Goeppert  Nico Goldscheider 《Hydrogeology Journal》2018,26(7):2241-2255

Karst aquifers are characterized by a high degree of hydrologic variability and spatial heterogeneity of transport parameters. Tracer tests allow the quantification of these parameters, but conventional point-to-point experiments fail to capture spatiotemporal variations of flow and transport. The goal of this study was to elucidate the spatial distribution of transport parameters in a karst conduit system at different flow conditions. Therefore, six tracer tests were conducted in an active and accessible cave system in Vietnam during dry and wet seasons. Injections and monitoring were done at five sites along the flow system: a swallow hole, two sites inside the cave, and two springs draining the system. Breakthrough curves (BTCs) were modeled with CXTFIT software using the one-dimensional advection-dispersion model and the two-region nonequilibrium model. In order to obtain transport parameters in the individual sections of the system, a multi-pulse injection approach was used, which was realized by using the BTCs from one section as input functions for the next section. Major findings include: (1) In the entire system, mean flow velocities increase from 183 to 1,043 m/h with increasing discharge, while (2) the proportion of immobile fluid regions decrease; (3) the lowest dispersivity was found at intermediate discharge; (4) in the individual cave sections, flow velocities decrease along the flow direction, related to decreasing gradients, while (5) dispersivity is highest in the middle section of the cave. The obtained results provide a valuable basis for the development of an adapted water management strategy for a projected water-supply system.  相似文献   

Hybrid Gaussian-cubic radial basis functions for scattered data interpolation     

Pankaj K. Mishra  Sankar K. Nath  Mrinal K. Sen  Gregory E. Fasshauer 《Computational Geosciences》2018,22(5):1203-1218

Scattered data interpolation schemes using kriging and radial basis functions (RBFs) have the advantage of being meshless and dimensional independent; however, for the datasets having insufficient observations, RBFs have the advantage over geostatistical methods as the latter requires variogram study and statistical expertise. Moreover, RBFs can be used for scattered data interpolation with very good convergence, which makes them desirable for shape function interpolation in meshless methods for numerical solution of partial differential equations. For interpolation of large datasets, however, RBFs in their usual form, lead to solving an ill-conditioned system of equations, for which, a small error in the data can cause a significantly large error in the interpolated solution. In order to reduce this limitation, we propose a hybrid kernel by using the conventional Gaussian and a shape parameter independent cubic kernel. Global particle swarm optimization method has been used to analyze the optimal values of the shape parameter as well as the weight coefficients controlling the Gaussian and the cubic part in the hybridization. Through a series of numerical tests, we demonstrate that such hybridization stabilizes the interpolation scheme by yielding a far superior implementation compared to those obtained by using only the Gaussian or cubic kernels. The proposed kernel maintains the accuracy and stability at small shape parameter as well as relatively large degrees of freedom, which exhibit its potential for scattered data interpolation and intrigues its application in global as well as local meshless methods for numerical solution of PDEs.  相似文献